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Abstract 

It is shown that, for stationary isotropic turbulence, Taylor's well known dissipa- 
tion surrogate Du' s /L can be derived directly from the Karman-Howarth equation 
and is in fact a surrogate for inertial transfer, which becomes equal to the dissi- 
pation, as the Reynolds number tends to infinity. The expression found for the 
dissipation rate e is 



A 3 u 
L 



/3 r 



1 A 2 



where the coefficients A 2 and A% depend on the second- and third-order structure 
functions respectively and Rl = u'L/u is the Reynolds number based on the integral 
length scale L. Further, consideration of the spectral energy transfer processes 
shows that the dissipation rate is entirely determined by the energy injection rate. 
The role of the viscosity is merely to determine the way in which the turbulent 
system adapts to increasing injection rates, both by increasing the volume of k- 
space and by changing the shape of the energy spectrum. The quantity u ri /L is 
a measure of the inertial transfer rate and only becomes equal to the dissipation 
when the Reynolds number is large enough to permit scale-invariance. It is noted 
that similar conclusions hold for shear flows, such as Poiseueille flow, both laminar 
and turbulent, where the analogue of energy injection is the rate of doing work 
expressed in terms of the pressure gradient and the mean velocity. 
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1 Introduction 



The idea that the rate of kinetic energy dissipation per unit mass of fluid turbulence has 
something anomalous about it, is widespread: see PQ-[I] and references therein. It arises 
from a result established in 1935 by Taylor [5], who showed, essentially by dimensional 
arguments, that the dissipation rate could be written as proportional to u' 3 /l where u! is 
the rms velocity of the fluid and I is some length scale of the system. Later, Batchelor [6] 
discussed this idea further in terms of the decay of kinetic energy (of turbulent motion) 
and offered two interpretations of it. The first of these was to see it as the decay of an 
amount of energy u' 2 in a time l/u'. His second interpretation was to regard it as the 
effect of an eddy viscosity ul acting on a shear of order u/l to 'produce a "dissipation" 
(sic) of energy from the energy- containing eddies to smaller eddies'. It seems likely from 
the context, and his use of quotation marks, that Batchelor saw this kind of expression 
as an approximation to inertial transfer which would be equal to the actual dissipation, 
under conditions of local statistical equilibrium. He also made use of the integral length 
scale, while noting that it was not as directly representative of the energy containing 
eddies as it might be. 

Essentially this is a pragmatic choice and in recent years the integral length scale L 
has been generally employed. Taylor's expression is normally written as 

e = Du' 3 /L, (1) 

where D tends to a constant^ with increasing Reynolds number [B]-[I2], and is widely 
used as a surrogate expression for the dissipation. 

Taylor's analysis suggested that D would depend on the the geometrical nature of the 
boundary conditions; while Batchelor pointed out that it might depend on the time of 
decay, the initial conditions of the turbulence, and the choice made for /. This view has 
received recent theoretical support [13], [13]. 

Now let us turn our attention to the so-called dissipation anomaly. Given that dissi- 
pation is due to viscosity and that the dissipation rate is formally defined in terms of the 
coefficient of viscosity, thus: 

the perceived anomaly takes one of two related forms: 

1. The fact that the dissipation rate, as given by ([T]) and verified by experiment, is 
found to be independent of the fluid viscosity [7] -[12], providing that the Reynolds 
number is sufficiently large. 

2. The existence of finite dissipation in the limit of vanishingly small viscosity (or the 
limit of infinite Reynolds number) pE]-[3]. 

Noting that ([T]) has only been derived from dimensional considerations, and analysed 
qualitatively in terms of the kinetic energy u' 2 and the eddy turnover time L/u', there 
is a need for a more theoretical analysis. From our present point of view, the recent 
work of Doering and Foias [15] is of interest. For the case of forced turbulence, they have 

1 Some workers in the field use C e for this constant. 



2 



established both upper and lower bounds on the dissipation rate. We shall return to this 
work at appropriate points in our own analysis. 

In this paper we mainly confine our attention to stationary, isotropic turbulence. We 
examine the exact relationships expressing conservation of energy, first in real space, and 
then in wavenumber space. We show that Taylor's expression follows quite naturally 
from the Karman-Howarth equation, as the Reynolds number increases. In the process, 
we find that it is a surrogate for the inertial transfer, and not for the dissipation rate, 
as such. Then, by considering the equivalent fc-space relation, we verify that the rate of 
doing work by body forces is the controlling quantity which determines the dissipation 
rate in this kind of turbulence; and by, fairly obvious extension, in all fluid flows. 



2 Taylor's surrogate and the Karman-Howarth equa- 
tion 

The Karman-Howarth equation may be written in terms of structure functions as [16J: 

_2 ids 1 = ±d_ (r , s ) _ ( dSA 

3 2 dt 6r 4 dr V V r 4 dr \ dr J ' 1 ' 

where the structure function of order n is given by 

S n = ((u(x + r)-u(x)) n ). (4) 

In this section we first consider the stationary case; then the relationship of the work 
we present here to other recent work on the subject; and finally the extension of our 
analysis to freely decaying turbulence. 



2.1 Stationary turbulence 

For the case of stationary turbulence, we may set the time-derivative equal to zero, and 
re-arrange the Karman-Howarth equation to obtain an expression for the dissipation, 
thus: 

1 d 4 3v d ( 4 dS 2 \ , , 

£ ~ ^d~r {rbz) + Wd~rX^r~)- {b) 

Now make the change of variables, 

S n = V n f n (x) J with x = T -, (6) 

where V is some constant velocity scale and b is any length scale. We should note that 
this step involves no approximations or non-trivial assumptions. It merely introduces the 
/„ as dimensionless forms of the structure functions. As they are dimensionless, their 
dependence on r must be scaled by some length, here denoted by b. In Section 2.3 we 
will discuss the introduction of self-similar and similarity solutions when we extend the 
present analysis to time-dependent (i.e. freely decaying) turbulence. 
With these substitutions, equation (0) becomes: 

A 3 V S A 2 vV 
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where the coefficients A 3 and A 2 are given by 



A, 



J_d_ r 

4a; 4 9a; 



x A h{x) 



and 



A, 



2x 4 <9x 



x 



dx 



(9) 



Taking out common factors, it is readily seen that the expression for the dissipation 
becomes 

£= — i 1+ ivd' (10) 

where the Reynolds number is given by Rb = Vb/v. We note that this equation is still 
just the Karman-Howarth equation: no approximation has been made. 

Now we may identify the prefactor on the right hand side of ( flOl) as being the same 
as the Taylor dissipation surrogate (0Q); providing we make the choices V = u', the root- 
mean-square velocity, and b = L, the integral scale. Then we obtain the general relation 



A 3 u 



/3 



Rl A 3 . 



but now the Reynolds number is given by Rl = u'L/v. Note that the f n are determined 
by this choice and hence also the coefficients A 3 and A 2 . 

Clearly, as the Reynolds number goes to infinity, this expression reduces to Taylor's 
surrogate for the dissipation §T§, provided that A 3 becomes a constant. However, as it 
represents the inertial-transfer term, we should really describe it as Taylor's surrogate for 
inertial transfer. 



2.2 Comparison with other work 



It is also of interest to compare equation ([7]), which is an intermediate stage in our 
calculation, to the result for an upper bound on the dissipation as given by Doering and 
Foias |15j . This is featured in their abstract as 

u' 2 u' 3 
e<c 1 v— + c 2 — , 

and corresponds to their equation (40). Here the coefficients c\ and c 2 depend on the 
shape of the forcing function, while I is its longest length-scale. Obviously this is quite 
different from our own result, where the corresponding parameters depend on the fluid 
turbulence and not on the forcing, and we have an equality, rather than an inequality. 
Nevertheless, the general similarity of the two results is worthy of notice. 

After some further manipulations, these authors introduce a function (3 (the same as 
our D in (0Q)) such that 

where the symbols are all as in [15] . They then use the identity Re = (3R\ to substitute 
for Re, solve the resulting quadratic equation, and obtain 



1 + + 



4a 



b 2 R\ 



(13) 
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in terms of the Taylor- Reynolds number R\. 

An interesting development is that Donzis, Sreenivasan and Yeung pj] have taken 
this upper bound as an equality, which they write as 

P = A(l + y/l + (B/R X )^ . (14) 

They fit this curve to results obtained from a numerical simulation and obtain an im- 
pressively close fit, with A ~ 0.2 and B ~ 92, leading to an asymptotic value of (3 = 0.4. 

Evidently our equation (iTTj) . being readily reduced to the same form as (TT2"j) . may be 
further reduced to the same form as (1141) . with 

A = A 3 /2 and B = 2A l/2 /A 3 . (15) 

Thus the implication is that it also agrees well with the results from simulation. We are 
currently doing our own numerical calculations; but, in the meantime, it is reassuring to 
know this. 



2.3 Extension to freely decaying turbulence 

We now consider whether our results can also be applied to freely decaying isotropic 
turbulence. The neglect of the time-derivative term in ([3]) is usual, even for freely decaying 
turbulence, provided that the Reynolds number is large enough and one restricts attention 
to the inertial range. For instance, this step is required in order to derive the well-known 
'4/5' law and is known as local stationarity. 

Let us consider the effect on © of choosing V and b to be the Kolmogorov velocity 
and length scales, as given by: 

v = (ue) 1/4 and rj = (is 3 /e) 1/4 : (16) 

see equation (21.4) in [T7j or equations (2.131), (2.132) in [18]. For n = 2, substitute (TIE]) 
into (J6j for V and b, respectively. We obtain 

S 2 (r) = p^e^Ur/rf). (17) 

If we are to make any further progress, then we must assume a (specific) self-similar form 
for f2- That is, we assume: 

Mr/v) = ^Mr), (18) 
from which the second-order structure function becomes: 

S 2 (r)=e 2 / 3 / 2 (r); / 2 (r)=r 2 / 3 . (19) 
Similarly, we can show that the third-order structure function is of the form 

S 3 (r) = e/ 3 (r); Mr) = r. (20) 

It should be noted that the forms in ( flBT) omit constants of order unity and this is reflected 
in these results. Also, we do not consider the question of non-canonical exponents or 
intermittency corrections as these matters have been addressed elsewhere [T9"] . 

Evidently our analysis is consistent with the K41 picture and should apply for the case 
of local stationarity. Let us now consider a more general situation, where we assume that 
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the decaying turbulence is self-preserving. Now we replace (JSJ) by the time-dependent 
form: 

S n (r,t) = u' n (t e )f n (x), (21) 

where 

x = r/L(t). (22) 

Note that our choice of an initial time for a similarity solution requires some care. We 
have taken it to be some t = t e , when the turbulence is said to have evolved from 
arbitrary initial conditions. If we assume that the turbulence is self-preserving after 
t = t e , then there is no explicit time- dependence in the dimensionless structure functions, 
and equation ffTTj) applies in the present case as well. 

In practice self-preservation is only likely to be approximately correct, and we can 
take account of the residual time-dependence as a perturbation of the similarity solution, 
writing equation ( 12~TT) as 

S n (r,t) = u' n (t e )f n (x,r), (23) 

where now 

r = t/T; T = L(t)/u'(t e ). (24) 

In these circumstances it is easily shown that equation (Hip can be generalised to the 
form: 

(A 3 - B 2 )u' 3 L 1 A, 



L 

where the new coefficient B 2 is given by 



1 + 



Rl ^3 — B 2 



(25) 



for small values of the time derivative. 

In all, this suggests that equation (flTl) . although derived for stationary turbulence, 
may apply quite well to decaying turbulence. 



3 The spectral picture 

In order to examine these ideas further, we study the spectral energy balance in wavenum- 
ber space, and consider the effect of increasing the Reynolds number, ultimately taking 
the limit of infinite Reynolds numbers. The three energy transport processes, viz., injec- 
tion, dissipation and inertial transfer will now be considered in turn. 
We begin with the energy balance equation, in its well known form, 

^+2ujAE(k,t) = T(k,t) + W(k), (27) 

where the energy transfer spectrum T(k,t) is given by: 

T(k,t) =2nk 2 M a/3 ^k) J d 3 j {C^ a (j,k-j,-k,t)-C^ a (-j,-k + j,k,t)}, (28) 

and C aj a 7 (k,j,k — j) is the three- velocity correlation. We may write T(k,t) as 

T(k,t)= / S(k,j;t)dj, (29) 
Jo 
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where S depends on the triple moment: its form can be deduced from (j28p . It can be 
shown that S is antisymmetric under the interchange k ^ j: 

S(k,j;t) = -S(j,k;t). (30) 

Hence _ ™ ™ 

POO POO POO 

/ T{k,t)dk = dk dj S(k,j;t) = 0, (31) 
Jo Jo Jo 

is an exact symmetry which expresses conservation of energy. 

In order to study the stationary case, we have added an input spectrum W(k) (if one 

wishes, this can be related to the covariance of the random stirring forces [IE])- We also 

introduce Ew as the rate at which the stirring forces do work on the turbulent fluid: 

POO 

e w = / W{k)dk. (32) 



The dissipation rate Ed in fc-space is defined by Ed = —dE/dt for freely decaying 
turbulence. As is well known, we can obtain an expression which is also valid for the 
stationary case, in the usual way, by temporarily setting W(k) = 0, integrating ( 1271) over 
wavenumber, and rearranging, such that the energy balance becomes: 

dE r°° 

e d = — -=/ 2uk 2 E(k,t)dk, (33) 
dt Jo 

where we have also invoked equation ( 13T1) . The region in fc-space where the dissipation 
mainly occurs is characterised by the Kolmogorov dissipation wavenumber: 

k D = (£D/f 3 ) 1/4 - (34) 

Then, restoring the injection spectrum, for the stationary case we have dE(k,t)/dt = 
0, and the energy balance becomes: 

T(k) + W(k) - 2vk 2 E(k) = 0. (35) 

Integrating both sides with respect to wavenumber, we have: 

/■OO /"OO 

/ W(k)dk - / 2vk 2 E(k)dk = 0; or: e w = e D , (36) 
Jo Jo 

as a rigorous consequence of stationarity. In other words, for a stationary flow, the 
dissipation is governed by the rate at which we do work on the fluid in order to produce 
a required fluid motion. We shall enlarge on this point presently. 



3.1 Inertial transfer and scale- invar iance 

We now consider inertial transfer of energy in fc-space. The energy flux is introduced 
when we integrate each term of (1271) with respect to wavenumber, from zero up to some 
wavenumber k. Reverting for the moment to the general (i.e. non-stationary) case, we 
obtain: 

d t K r°° c K r K 

— dkE(k,t) = - dk dj S(k,j;t) - 2u dkk 2 E(k,t), (37) 
dt Jo Jk Jo Jo 

where we have used the antisymmetry of S, and made some rearrangements. In this 

form the effect of the transfer term is readily interpreted as the net flux of energy from 
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wavenumbers less than k to those greater than k, at any time t [BJ. Denoting this flux by 
n(/t); and, in order to avoid the ambiguity associated with the scale-invariance paradox, 
making an exact decomposition of the transfer spectrum into filtered-partitioned forms 
T+-(k\K) and T~+(%) [20], we have 

/"CO rK 

n(«)=/ dkT + -(k\n) = - dkT- + (k\n), (38) 

where we have now assumed stationarity and dropped the time dependence. (Note that 
the decomposition is completed by T {k\n) and T ++ (A;|k), which are separately conser- 
vative on the intervals [0,k] and [k, oo], respectively [20j . ) The maximum value of the 
energy flux is U max (K), where T + ~(k\n) = T~ + (k\K) = 0. 

When setting up stationary, isotropic turbulence by means of an arbitrary choice of 
stirring forces, it is usual to try to reproduce the characteristic features of the classic 
turbulent shear flows. In terms of the energy spectrum, these may be seen as the energy- 
containing range (low wavenumbers), the inertial range (intermediate wavenumbers) and 
the dissipation range (large wavenumbers). It has been known since the late 1930s [6 J 
that the energy- containing and dissipation ranges become more widely separated as the 
Reynolds number is increased; with the Kolmogorov dissipation wavenumber ko, as given 
by (|34|) . providing a reliable measure of this process. In practice, what this means is that 
we should choose the injection spectrum W(k) to have a somewhat peaked form at low 
values of wavenumber, when compared to kr>. 

We may formalize this picture as follows. At sufficiently large R\, the energy- 
containing and dissipation ranges become separated by the inertial range of wavenumbers, 
thus: 

hot < « < hop, (39) 

where k now stands for any wavenumber in the inertial range. In this case, the injection 
and dissipation spectra satisfy approximate relationships as follows: 

/•hot roc 

/ dkW{k) ~ e w ; and / dk 2uk 2 E(k) ~ e D . (40) 

JO J ktop 

In this range of wavenumbers the maximum energy flux should be approximately constant 
and we will find it helpful to introduce a specific symbol for this quantity, thus: 

n ma x = £T- ( 41 ) 

For stationarity, we must have the overall energy balance: 

e w = e T = e D . (42) 

We note that these three different physical processes are normally denoted by the single 
symbol e, this being justified by their all being numerically equal. In our view, it is 
necessary to draw a distinction between them in order to avoid confusion. 

It should now be apparent that Taylor's expression, as given by (JTJ, should be written 
as Et = Du' 3 /L; which also explains the observed dependence on Reynolds number (see 
Fig. 1 in [8J). In other words, ([T]) becomes equal to the dissipation when (|42p is satisfied. 
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3.2 The limit of infinite Reynolds numbers 

The separation of the energy-containing and dissipation ranges by a scale-invariant in- 
ertial range has an interesting consequence. It allows us to obtain separate low-fc and 
high-fc balance equations; by first integrating from zero up to kb ot , and then from infinity 
down to k top . First, we have 

/ dk djS(k,j)+ W(k)dk = 0. (43) 
Jo Jk bot Jo 

That is, energy supplied directly by the input term to modes with k < kbot is transferred 
by the nonlinearity to modes with j > k^ot- Thus T(k) behaves like a dissipation and 
absorbs energy. Second, 

roo rktop roo 

/ dk djS{k,j)- 2uk 2 E(k)dk = 0. (44) 

That is, the nonlinearity transfers energy from modes with j < kt op to modes with 
k > k top , where it is dissipated into heat. In this range of wavenumbers T(k) behaves like 
a source and emits energy which is then dissipated by viscosity. 

Now we reach the crux of our argument. We examine the limit of infinite Reynolds 
number, as follows. We keep the injection spectrum W{k) fixed (and hence also By/ = 
£t = £d are held constant), and allow the viscosity to tend to zero. As a result (originally 
pointed out by Batchelor in 1953 [6]), the sink of energy is displaced to k = oo, in the 
limit of infinite Reynolds number. We can deduce that this is so, either from the form 
of the Kolmogorov dissipation wavenumber or the form of the dissipation term; or even 
from a consideration of the local Reynolds number for mode k [S]. 

However, we cannot emphasise too strongly that this has nothing to do with the Euler 
equation. The Euler equation can indeed be obtained from the Navier-Stokes equation 
by setting the viscosity equal to zero. The result is the equation of motion for an inviscid, 
zero-dissipation fluid. Here, in contrast, we take the limit of zero viscosity, while keeping 
the dissipation constant. 

This operation was later formalized by Edwards [2T] (or see Section 6.2.7 in [IE]), who 
argued that the dissipation at infinity could be represented by a Dirac delta function, 
thus: 

roo roo 

/ dklim2uk 2 E(k) = dke D 5{k - oo) = e D . (45) 

Jo ^—o e D =constant Jo 

Similarly, with an infinite extent of wavenumber space, any injection spectrum of finite 
extent can be scaled back to the origin, so that we have also: 

roo roo 

/ dk lim W(k) = dke w 8{k)=e w . (46) 

Jo 1 e w =constant -^o 

Thus, in the limit of infinite Reynolds numbers, equation (|3"5"]) for the spectral energy 
balance may be written as: 

- T(k) — ew S(k) — e D 8(k — oo). (47) 



Note that these forms satisfy all the relevant relationships given above as equations ( 1391) 
- (jUJ) and, although this may seem a rather extreme procedure, it is in fact nothing more 
than a different mathematical representation of scale invariance in the inertial range. 
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The general forms of our arguments given here in spectral space are not unlike those 
of Davidson (see Section 3.2.2 in |22j), which are presented for physical or scale space. 
A particular point of interest is his example of a sudden expansion in a pipe, where the 
'head loss' is determined by continuity and the momentum theorem (plus, it must be 
said, an approximation based on the existence of stagnation points) but not apparently 
on the fluid viscosity. 

Also, it is of interest to note that our views expressed here seem to be quite close to 
those of Doering and Foias [IS], despite the very different form of the analysis. Essen- 
tially both treatments give a dominant role to the production of the turbulence and, for 
instance, our equation fl36l) is essentially just equation (17) of reference [15]. 

4 Conclusion 

From this simple analysis, it is apparent that there are really no grounds for viewing 
the lack of dependence of the dissipation on the viscosity as an anomaly. Conservation 
of energy ensures that the dissipation is just equal to the injection rate. The role of 
the viscosity is to control the energy occupation of wavenumber space, along with the 
shape of the spectrum. Noting that Taylor's expression for the dissipation, as given 
by equation (OQ), is really just the inertial transfer rate, we can appreciate that it will 
depend on Reynolds number, until the latter is large enough for the energy-containing 
and dissipation ranges of wavenumber to be adequately separated, to the point where 
inertial transfer is well-defined, and equal to the actual dissipation rate. Also, for the 
limit of infinite Reynolds number, it can be seen that the symmetry-breaking effect of 
the viscosity is still present in the form of the delta functions, as in equation (31). 

We can also attempt to say something more general about shear flows. If we consider 
classical shear flows, then laminar Poiseuille flow presents one of the few cases where 
the Navier-Stokes equation may be solved exactly. This is because the nonlinear term 
vanishes identically, and the problem resolves itself into a balance between viscous and 
pressure forces. Here the dissipation rate can be worked out exactly, and is found to be 
independent of viscosity. It is, as it must be, equal to the injection rate; which can also 
be worked out exactly, in terms of the rate at which the pressure force moves its point of 
application. In this case, the gradient of the velocity depends inversely on the viscosity, 
and becomes steeper in order to maintain the dissipation rate as the viscosity is decreased 
at constant axial pressure gradient. This is, of course, analogous to the behaviour of the 
Kolmogorov length scale under increasing Reynolds number at constant injection rate. 

Extension of the analysis to turbulent Poiseuille flow requires a little more work, but 
those nonlinear terms which do not vanish identically, do vanish when integrated over 
the system volume; and the analysis for the laminar case applies here as well. 
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